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Abstract. Reliability of a system is considered where the components' ran- 
dom lifetimes may be dependent. The structure of the system is described by 
an associated "lattice polynomial" function. Based on that descriptor, gen- 
eral framework formulas are developed and used to obtain direct results for 
the cases where a) the lifetimes are "Bayes-dependent" , that is, their inter- 
dependence is due to external factors (in particular, where the factor is the 
"preliminary phase" duration) and b) where the lifetimes' dependence is im- 
plied by upper or lower bounds on lifetimes of components in some subsets 
of the system. (The bounds may be imposed externally based, say, on the 
connections environment.) Several special cases are investigated in detail. 



1. Introduction 

A semicoherent system consisting of nonrepairable components with random life- 
times can be associated with a Boolean function called the structure function of the 
system. It expresses the "on" indicator of the system through the "on" indicators 
of the components based on the logic of connections. In turn, the structure function 
can be extended to a lattice polynomial function called the life function of the sys- 
tem, which expresses the system lifetime in terms of the component lifetimes using 
minimum and maximum in place of conjunction and disjunction, respectively. 

In this paper we use the results on lattice polynomial (l.p.) functions from 
a series of our recent articles to obtain formulas for system reliability where the 
components' lifetimes are dependent random variables. 

In the current literature, including such prominent sources as [21 [21 [TOl [17] , that 
case is explored only under the additional assumption of exchangeable components 
(Barlow and Proschan [3]) which presupposes that the marginal distributions 
of components' lifetimes are identical. Meanwhile, in many realistic systems that 
assumption cannot be applied. 

Using l.p. functions we derive a general formula for system reliability that yields 
specific results for some important special cases of dependence among components. 

One is the case where the dependence can be ascribed to the simultaneous in- 
fluence on all components of some random factors (such as, say, the system's envi- 
ronment physical parameters: pressure, temperature, moisture level, etc). For any 
given set of values of those factors components' lifetimes are conditionally indepen- 
dent (but their probability distributions depend on the factors' values). We will 
refer to this case as "Bayes-dependence" . 
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In particular, of special interest is a model where there is a preliminary period 
( "pre-phase" ) of random duration and all components are guaranteed to survive 
it. Once it is over, the components' individual "after-phases" are independent and 
their probability distributions depend on the "pre-phase" duration. 

Another pattern of dependence arises when there are collective upper or lower 
bounds on lifetimes of certain subsets of components imposed by external conditions 
such as physical properties of the assembly or its environment. Say, when some 
components of the system are connected to the same power source, the lifetime of 
the source (or the fuse) becomes a collective upper bound for those components' 
lifetimes. Or, when a part of the system is backed up by a reserve standby device 
(too expensive, perhaps, for a regular duty) that instantly picks up duties of any 
failed component in that subset, the device's lifetime becomes a lower bound for the 
components' lifetimes in that subset. A "senior" air traffic controller could serve 
as an example here. 

To our knowledge, models with such patterns of components' interdependence 
have not been analyzed in the literature. Using l.p. description of a system in this 
paper presents a natural framework for such analysis bcicause each component's 
service duration can be easily and naturally represented using an l.p. function in- 
volving the component's own lifetime and bounding random variables. 

Moreover, the case where the collective bounds arc constant can be analyzed 
through so-called weighted lattice polynomial (w.l.p.) functions. 

This paper is organized as follows. In Section 2 we introduce the l.p. function of 
a semicoherent system and use it to present a number of reliability formulas for a 
system with generally dependent components' lifetimes. 

In Section 3 we describe the case of "Bayes-dependence" of components' lifetimes 
and present exact formulas for the system reliability. In turn, those formulas make 
it possible to provide exact formulation of reliability parameters such as the mean 
time-to-failure of the system. 

In Section 4 we analyze systems with "pre-phase" dependent components. In 
particular, we present closed-form results for the case where the "after-phase" du- 
rations of components are assumed to have exponential distributions with failure 
rates depending on the "pre-phase" duration. We consider this model as giving a 
system analyst an alternative to using Weibull distributions, more complex both 
analytically and statistically. 

In Section 5 we consider systems with lower and/or upper bounds on service 
duration of (subsets) of components. There, we show how the system's reliability 
can be computed as reliability of another, "augmented" system with added compo- 
nents and connections but without subset bounds. Based on that representation, 
when the original components intrinsic lifetimes are independent, the augmented 
system's reliability is given by formulas of Sections 2 and 3. 

For any numbers a. ft G K = [— cx3, cxd] and any subset A C [n] = {1, .... n}, let 
e^'^ denote the characteristic vector of A in {q,/3}", that is, the n-tuple whose ith 
coordinate is ft, if i € A, and a, otherwise. 
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2. Lattice polynomial function and reliability of a system 

Consider a system consisting of n components that are interconnected. The state 
of a component i g [n] can be represented by a Boolean variable Xi defined as 

{1, if component i is functioning, 
0, if component i is in a failed state. 

For convenience, we also introduce the state vector x — (xi, . . . , x„). 

As is common in the literature, the state of the system is described from the com- 
ponent states through a Boolean function (p: {0, 1}" — > {0, 1}, called the structure 
function of the system and defined as 

I 1, if the system is functioning, 
0(x) — < 

I 0, if the system is in a failed state. 

We shall assume throughout that the structure function is nondecreasing (the 
system is then said to be semicoherent) and nonconstant, this latter condition 
ensuring that 0(0) = and 0(1) = 1, where = (0, . . . , 0) and 1 = (1, . . . , 1). 

As a Boolean function, the structure function can also be regarded as a set 
function v. 2["1 — > {0, 1}. The correspondence is straightforward: We have v{A) — 
0(e^ ) for all A C [n] and 

(1) 0(x)= j2 «(^)n^^ n 

A(Z[n] ieA ii£[n]\A 

We shall henceforth make this identification and often write 0«(x) instead of 0(x). 
Clearly, the structure function 0^ is nondecreasing and nonconstant if and only if 
its underlying set function v is nondecreasing and nonconstant. 

Another concept that we shall often use in this paper is the dual of the set 
function v, that is, the set function v* : 2["1 — ;> {0, 1} defined by 

v*{A) = l-vi[n]\A). 

For any event E, let Ind(£') represent the indicator random variable that gives 1 
if E occurs and otherwise. For any i G [n], we denote by Ti the random time-to- 
failure of component i and we denote by Xi{t) — Ind(ri > t) the random state at 
time t ^ of component i. For simplicity, we introduce the random time-to- failure 
vector T = (Ti, . . . , T„) and the random state vector X(t) = {Xi (i), . . . , X„(t)) at 
time t ^ 0. We also denote by Ts the random time-to-failure of the system and by 
Xs{t) = Ind(T5' > t) the random state at time i > of the system. 

The structure function clearly induces a functional relationship between the 
variables Ti, . . . , T„ and the variable Ts- As we will see in Theorem [21 Ts is always 
an l.p. function of the variables Ti, . . . , T„. Just as for the structure function, this 
l.p. function provides a complete description of the structure of the system. 

Let us first recall the concept of l.p. function of real variables; see for instance 
Birkhoff [4] and Gratzer [8]. Let L C M. denote a totally ordered bounded lattice 
whose lattice operations A and V are respectively the minimum and maximum 
operations. Denote also by a and b the bottom and top elements of L. We assume 
that a ^ b. 

Definition 1. The class of lattice polynomial (l.p.) functions from L" to L is 
defined as follows: 
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(i) For any k £ [n], the projection t = (ti, . . . , t„) i-^ tk is an Lp. function from 
L" to L. 

(ii) If p and q are Lp. functions from L" to L, then p A q and p \/ q are Lp. 
functions from L" to L. 

(iii) Every Lp. function from L" to L is constructed by finitely many appUcations 
of the rules (i) and (ii). 

Clearly, any Lp. function p: L" — >■ L is nondecreasing and nonconstant. Fur- 
thermore, it was proved (see for instance [5) that such a function can be expressed 
in disjunctive and conjunctive normal forms, that is, there always exist noncon- 
stant set functions w"^: 2["1 — > {a,b} and tu^: 2["1 with w'^{0) — a and 
w'^{0) = b, such that 

(2) Pit)= V A*»= A V*- 

AC[n] i£A AC[n] ieA 

w<'{A)=b w'={A)=a 

Clearly, the set functions w'^ and w'^ that disjunctively and conjunctively define 
the Lp. function p in ([5]) are not unique. However, it can be shown [TT] that, from 
among all the possible set functions that disjunctively define p, only one is non- 
decreasing. Similarly, from among all the possible set functions that conjunctively 
define p, only one is nonincreasing. These special set functions are given by 

w''{A)=p{ef) and ^=(A) = p(e['4^). 

The Lp. function disjunctively defined by a given nondecreasing set function w : 2["1 — 
{a, b} will henceforth be denoted Pw We then have 

put)= V A^*= A V^- 

AQln] ieA AC[n] ieA 

w{A)=b w*{A)=b 

where w* is the dual of w, defined as w* =70 (7^^ o w)* , the function 7 : {0, 1} — )■ 
{a, 6} being a simple transformation defined by 7(0) — a and 7(1) — b. 

Remark 1. From any nonconstant and nondecreasing set function w: 2^'"^ — > {a, b}, 
define the set function : 2["1 {a, b} as 



b, if w{A) = b and w{B) = a for aU B <^ A, 
a, otherwise. 



The disjunctive and conjunctive representations of the Lp. function p^, having a 
minimal number of terms are given by 

pv.it) ^ V f\*^= A V*^ 

AC[n] ieA AC[n] ieA 

u„(yl)=b u„.(A)={) 

(see Proposition 8 in [11]) and are exactly those minimal paths and cuts represen- 
tations of Pw (as given, say in Theorem 3.5 of |2j), namely 

r s 

p-(t) = V A = A V 

3=itePj 3=iieKj 

where Pi , . . . , are the minimal path sets and Ki , . . . , Kg are the minimal cut 
sets. 
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The following theorem points out the one-to-one correspondence between the 
structure function and the l.p. function that expresses Ts directly in terms of the 
variable Ti, . . . , r„. It is the latter fact that makes, in our opinion, the l.p. function 
a preferable system descriptor. 

As lifetimes are [0, cxD]-valued, we shall henceforth assume without loss of gen- 
erality that L = [0, oo], that is, a = and b — oo. Recall that the coproduct of n 
Boolean variables defined by 11^2;^ = 1 — 11^(1 — Xi). 

Theorem 2. Consider an n-component system whose structure function (p is non- 
decreasing and nonconstant. Then we have 

(3) Ts =PwiTi,...,Tn), 

where w = ^ o v. Conversely, any system fulfilling ^ for some l.p. function 
Pw- i" L has the nondecreasing and nonconstant structure function (jjy, where 
V — 7^""" o w. 

Proof. The proof mainly lies on the distributive property of the indicator function 
Ind(-) with respect to disjunction and conjunction: for any events E and E' , 

lndiE\JE') = lnd{E)\Jlnd{E'), 

lnd{EAE') = Ind(i;) Alnd(£;'). 

Thus, for any t ^ we have 

Ind(p^(T)>i) = Ind( V >t) = V /\lnd{T,>t) 

AC[n] ieA A(l[n] ieA 

v{A) = l v{A) = l 

= n n^w 

AC[n] ieA 
v(A) = l 

Hence, we have Ts — Pw{T) if and only if Xs{t) — 0„(X(i)) for all t ^ 0, which 
completes the proof. □ 

The reliability function of component i is defined, for any t ^ 0, by 

R,it) = Pr(r, > = Pr(X,(t) = 1) = E[X,it)], 

that is, the probability that component i does not fail in the time interval [0,t]. 
Similarly, for any t ^ 0, the system reliability function is 

Rs{t) - Pr(r5 >t)^ Pr(Xs(t) - 1) = nXs{t)l 

that is, the probability that the system does not fail in the time interval [0,t]. 

Based on representation ([T]) and its dual form, we present general formulas for 
the system reliability function in case of generally dependent variables Ti , . . . , T„ 
(first presented in [B]). 

Theorem 3. We have 

(4) Rsit) = J2 viA)PriXit)=eY), 

AC[n] 

(5) Rsit) = 1- 5] «*(A)Pr(X(0=ej'4j. 

Adln] 



6 



ALEXANDER DUKHOVNY AND JEAN-LUC MARICHAL 



Name 


<l>v (x) 


Primal form 


E «(^) n n (1 - 

AC[n] ieA ie[n]\A 


Dual form 


1- j: v*{a) n n(i-^.) 

AC[n] i£[n]\A iGA 


Primal Mobius form 


E rriyiA) n 

AC[n] i&A 


Dual Mobius form 


E rny*{A) U Xi 

AQln] ieA 


Disjunctive normal form 


u viA) n X. 

AQ[n] ieA 


Conjunctive normal form 


n v*iA) u Xi 

AC[n] i£A 



Table 1 . Various forms of the structure function 



Proof. By ([T]), we have 

(6) Rs{t) = E[Mm)] = 5] i'(A)E[n^.W U (l-Xdt)) 

AC[n] ieA ie[n]\A 

= ^ v{A) Pr(X(i) = e:;^^), 

AG[n] 

which proves Formula ([SJ can be proved similarly by using the dual form of 
4>v (i.e., the second expression in Tabled]). □ 

More reliability function formulas can be obtained based on other structure func- 
tion representations. Any Boolean function has a unique expression as a multilinear 
function in n variables, 

(7) 0^(x) = J2 m4A)Y[xi 

ACln] ieA 

(see for instance [9 ), where the set function rriy : 2["1 — >■ Z is the Mobius transform 
of V, defined by 

m„(A)= ^(-l)l-^l-l^^(B). 

BCA 

By using the dual set function v* we can easily derive further useful forms of 
the structure function. Table [T] summarizes the best known forms of the structure 
function. 

Remark 2. Since is a Boolean function, we can always replace in its expression 
each product 11 and coproduct 11 with the minimum A and the maximum V, re- 
spectively. Thus, Theorem [5] essentially states that is also an l.p. function that 
has just the same max-min form as but applied to binary arguments. More 
precisely, <pv is similar to pw in the sense that 7 o 0^ = p^, o (7, . . . , 7). 
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Consider the joint distribution function and the joint survival function, defined 
respectively as 

F(t) = Pr(ri ^ Vi e [n]) and i?(t) = Pr(T, > Vi G [ti]). 

By using the same argument as in the proof of Theorem [3l we obtain two further 
equivalent expressions of Rs (t) . 

Theorem 4. We have 

Rs{t) = m.(.A)R{eY) 



AC\' ^ 

]\A> 



Rs{t) = 1- ^ m,-(A)F(e*'°° 



ACl 



Proof. By d?]), we have 

Rs{t) = E[M^m = 51 "*"(^)E[n^*w 



A(l\n] ieA 



J2 my {A) Pr(T, > iVi e A) 

AQ[n] 

J2 my{A)Rie'/). 



AC\n] 



Similarly, using the dual Mobius form of (f)y (i.e., the fourth expression in Table [IJ, 
we have 



Rs{t)=E[MMt))] = E 

AQ[n] 


{A)e[1[xs)] 

ieA 


ACln 


m,.(A)E[[](l-X,(t)) 

ieA 


= 1-E 

AC[n 


m,,, {A) Pr(Ti ^ i Vi G A) 


= 1-E 

ACln 


my.{A)F{e\;°^^^), 


where we have used the fact that X^ACfn] "^f* 





□ 

Remark 3. Based on the minimal path/cut sets representations (see Remark [T]), 
one can find in Ij and further in, say, [S] and [13] linear representations of Rs{t) in 
terms of reliability functions "series" and/or "parallel" subsystems. In this paper 
we chose to use the representations of Theorem [3l among other reasons, in order to 
use the Mobius transform and obtain two further equivalent expressions of Rs(t). 

The mean time-to-failure of component i is defined as MTTF^ = E[Ti] and sim- 
ilarly the mean time-to-failure of the system is defined as MTTF5 = E[T5]. These 
expected values can be calculated by the following formulas (see for instance [T^) 

poo POO 

MTTFi / R,{t)dt and MTTFg = / Rs{t)dt. 
Jo Jo 
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It is noteworthy that Theorem |4] immediately provides concise expressions for 
the mean time-to-failure of the system, namely 



MTTFs = 
MTTFs = 




Theorem H] may suggest that the complete knowledge of the joint survival (or 
joint distribution) function is needed for the calculation of the system reliability 
function. Actually, as Theorem [3] shows, all the needed information is encoded in 
the distribution of the indicator vector X(t). In turn, the distribution of X(t) can 
be easily expressed (see [6l [7]) in terms of the joint distribution function of X(t), 
which is defined by 

Pr(X(0=e:ii)= 5](-l)I^H^I^^(er). 

BCA 

That is, the joint distribution function needs to be known only where its arguments 
are equal to i or oo. 

In the case when Ti, . . . ,r„ are independent, which implies that the indicator 
variables Xi{t), . . . ,Xn{t) are independent for all t ^ 0, from ^ we immediately 
obtain (see for instance jTS]): 

(8) i?s(i)= E ^(^)n^»w n (i-^^w)- 

AC[n] ieA ieln]\A 

By extending formally the structure function (j)^ to [0, 1]" by linear interpolation, 
we define the multilinear extension of 4>y (a concept introduced in game theory by 
Owen that is, the multilinear function (j>^ : [0, 1]" [0, 1] defined as 

(9) ^„(x)= E ^(^)n^^ n (i-^*)- 

AC[n] ieA ie[n]\A 

We then observe that each of the alternative expressions of introduced in 
Table [T] can be formally regarded as a function from [0,1]" to [0,1], which then 
identifies with the multilinear extension of 

Combining ([5]) and we retrieve the classical formula (see for instance [TBI 

Rs{t)^<f>,{Riit),...,Rn{t)) 

where the function 0^ is called the reliability polynomial. Thus, both Rs{t) and 
MTTFs can be expressed in different forms, according to the expressions of 0„ 
corresponding to Table [T] For instance, using the primal Mobius form of cf)^, we 
obtain 

AC[n] i£A 



MTTF, 




A similar expression was obtained for the general case in Section 5 of |14| . 
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Consider the special case when the reliability of every subset A C [n] depends 
only on the number |^| of components in A (which happens, for example, when the 
component lifetimes are exchangeable). That is, 

R{e°/) = R{e"x') whenever \A\ = \A'\, 

and similarly for F(e*^^^). 

Defining R{k,t) = R{e°/), F{k,t) = F(e[;^'j^^), where k = \Al and 

my{k) = ^ TO-u(^), 

A<Z[n\ : \A\=k 

from Theorem 2] we derive the following corollary^ 

Corollary 5. // the reliability of every subset depends only on the number of com- 
ponents in the subset, then 

n 

Rsit) - ^rny{k)R{k,t) 

k=l 

n 

Rs{t) = l-^m„.(fc)F(fc,i). 

fc=i 

3. Systems with "Bayes-dependent" component lifetimes 

Consider a system where components' functioning is influenced by certain (per- 
haps, random) factors that may be both internal and external to the system. Say, 
failure rates of individual components may be influenced by the system's current 
environment conditions, such as temperature, pressure, precipitation, etc. Or, the 
whole assembly is made of subsystems that have "central" units and their status 
directly affects the other units of the subsystem. 

One way to model that situation is to introduce a set of (random) factors 
[/i, . . . , Urn whose joint probability density function g{ui, . . . , Um) defined on a do- 
main Dg is known. Given a fixed set of the factor's values, the system's components' 
lifetimes are assumed independent with individual conditional cumulative distribu- 
tion functions Fi{t,ui, . . . ,Um) depending on the factors' values. In that case, the 
joint probability distribution function of the components' lifetimes acquires an "in- 
tegrated product" form: 




We will refer to such interdependence pattern as "Bayes-dependence" . 

Clearly, when values of factors Ui, . . . , Um are fixed, since Ti, . . . , are (con- 
ditionally) independent, the indicator variables Xi{t),. . . ,Xm{t) are too, so using 
([T0| in (IH) with fixed u the same way as in the derivation of ^ we then integrate 



^For the case of exchangeable component Ufctimcs, this corollary was actually obtained in |13| 
by using the concept of signature. 
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over all possible values of u and obtain a generalization of (jH} : 
Rsit) = [ g{u)Pr{Ts>t\u)du 

JDg 

= I g{vi) Y v{A)\{R,{t,M) J] (l-i?,(i,u))du, 

AC[n] iGA jG[ii]\A 

where i?i(t, u) = 1 — Fi{t, u). 

In turn, it leads to a generalization of the classical reliability formula 

Rs{t)^ I giu)J,{Ri{t,u),...,Rn{t,u))du. 

Once again, both Rs{t) and MTTFg can be expressed in different forms, ac- 
cording to the expressions of (j)^ corresponding to Table [T] In particular, it is very 
convenient for calculations to use the primal Mobius form of cf)^: 

(11) Rsit) = ^ m,(^) / g(u) J]i?,(i,u)du, 

AC[n] ^9 leA 

(12) MTTFs = Y '^^i^) I / \{Ri{t,\i)dtdvi. 

In some cases it is natural to think of the "intrinsic randomness" of a unit's life- 
time (expressed in the shape of its distribution) as specific to the unit itself, while 
the overall parameters of the lifetime (mean, variance, etc.) can be influenced by 
the external factors. One way to model that is to consider only the distribution 
parameters as functions of those factors (without changing the form of the distri- 
bution) . 

For example, the assumption of constant unit failure rate is widespread in the 
literature and often justified by data. At the same time, external factors, physical or 
societal, often influence values of the rates. In the framework of Bayes-dependence, 
one can model that situation by setting Ri{t, u) = g-->>i(u)t^ 

Computing the integrals in the formulas above, one obtains: 



Rs{t) = "'-(^) / e-^-(")*5(u)du, 

AC[n] ■'^a 

MTTF5 = V m,{A) f du, 
where Aa(u) = J^ieA ^li'^)- 

Remark 4. The above results actually require merely that the joint distribution 
function of the components' lifetimes should have the integrated product form of 
([TU)l only when the arguments of the distribution function are either equal to t or 
oo. 

4. Systems with pre-phase 

In practice and literature many cases emerge where the units' failure rates should 
be modeled as nonconstant, say, increasing with elapsed time. A popular way to 
address that issue is to present the rates as functions of the elapsed time using, 
for example, WeibuU distribution. However, so constructed models are, for obvious 
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reasons, much harder to analyze mathematically and provide reliable estimates 
of parameters statistically. Moreover, in the context of reliability systems it is 
impossible to ignore the influence of joint functioning and interaction on individual 
units' failure rates. 

The framework of Bayes-dependence yields a natural way to make a model where 
both the joint functioning and elapsed time are accounted for and the analytical 
convenience of exponential distribution can be retained. 

Let each component's lifetime Tj consist of a random variable ( "pre-phase" ) U 
common for all components followed by, once the pre-phase is over at some time 
u, the subsequent individual "decay" phase Yi. Thus, both the variability of the 
system's failure rate and the in-system interaction effect can be modeled through 
the pre-phase distribution, leaving the analyst enough room to model the residual 
lifetime of the system separately. 

We will denote by G(u), g{u), and E[U] the distribution and density functions 
and expectation of U . The conditional distribution function of each decay phase 
given that U = u will be denoted by Fi(y, u). 

This case, too, belongs in the class of Bayes-dependence, where the pre-phase 
duration U plays the role of an external factor and 

(13) R,{t,u) = l-F,{\t-u\ + ,u), 

where we denote = max(a;,0). 

Substituting (fTB]) in (fTT|) and using the identity J2AC[n] "^til^) = 1, we obtain 

poo 

^s{t) = ^ m^{A) g{u)Y[R,{t,u)du 

AC[n] " ieA 

/oo pt 
g{u)du+ ^ my{A) / g{u) Y\_ {'^ - Fi{t - u,u)) du 

A^lnl - AC[n] "^^ ieA 

= l-G{i)+ I ^ m,(A) Jl 

Adln] ieA 



Thus, we have 



where 



Rs{t) = I - G{t) + f R*s{t-u, 
Jo 



u) g{u) du, 



R*s{y,u)^ J2 m,{A)l[{l-F,{y,u)). 

AC[n] ieA 

From this result we obtain 

poo 

MTTFs = / Rs{t)dt 
Jo 

/•oo /'OO /'t 

= / {1-G{t))dt+ / R*g{t-u,u)g{u)dudt 
Jo Jo Jo 

poo poo 

= E[U]+ / / R*s{t-u,u)dtg{u)du 



OO poo 



E[U] + / / R*g{t,u)dtg{u)du 
Jo Jo 
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and hence 



where 



POO 

MTTFs = E[U] + / MTTFs(u) g{u) du, 
Jo 

1*00 POO 

MTTF5(u) = / R*s{t,u)dt^ ^ to„(A) / Y[ {l ~ F^{t,u)) dt. 



AC[n] iGA 

When all individual decay phases have exponential distribution with parameters 
Xi{u), that is Fi{y,u) = 1 — e"'*''^"^'', the formulas above reduce to 

Rsit) = 1-G(t)+ y m^iA) f e-^-^(")l*-"l+g(u)du, 



AC[n] 



MTTFs = E[U] + y m,{A) f 

Ann] 



Aa(u) 



where Xa{u) = J^ieA^ii'^)- 

The previous formulas provide analytically simple results when the pre-phase' 
distribution is modeled using a piecewise constant density function. For example, 
in the simplest case where the pre-phase is uniformly distributed in an interval [a, b] 
and decay failure rates are constant we have 

1, if t < a, 

+ > m^A) — — ^, if a ^ t =^ 6, 

Rs{t)^{^-^ AT{n] >^A{b-a) 

5. Systems with collective bounds 

The environment in which a system is installed typically imposes its constraints 
on the functioning of the system. A power source that feeds several components 
imposes an upper bound on their service duration - it becomes a maximum of the 
component's own lifetime and the source's one. 

There is extensive reliability literature on a model where upper bounds are im- 
posed by fatally disabling external shocks to individual components and subsystems. 
The shock emergence processes are assumed Poisson, independent of the system and 
independent for different components and subsystems. It has been shown that in 
the reliability terms the action of shocks on the system is equivalent to assuming 
a Marshall-Olkin type distribution for the component lifetimes. (One can see, e.g., 

m) ' 

In the context of insurance applications, a contract can be purchased with an 
on call emergency service that backs up services of some components. It provides 
a lower bound on their service duration, that is, component's service duration 
becomes a minimum of its intrinsic lifetime and the emergency service's one. 

Today, in a complex system component services may have a complex combination 
of upper and lower bounds imposed on them by the nature of the system's purpose 
and its environment. In general, the presence of collective bounds and interaction 
of components with them, all the maximums and minimums that emerge that way, 
can be represented as follows. 
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Denote by Qi, Q2, ■ • • , Qm bounding random variables and by - the intrin- 
sic lifetime of the ith component (if left alone), the component's service duration 
emerging as a result of interaction of its intrinsic lifetime and relevant bounding 
factors can be expressed using an l.p. function qi 

T, = q,{T^,Qu...,Q^). 

Using that in Theorem [5] and p as the l.p. function of the system, the overall 
system's time-to-failure can be written as 

(14) Ts = p{qi{Tl Qi, . . . , Q„), . . . , g„(r°, Qi, . . . , Q™)). 

Remark 5. When all bounding variables are constant equation presents T5 as 
a so-called "weighted" lattice polynomial function of intrinsic lifetimes; see 

A composition of l.p. functions, it follows immediately from the definition, is an 
l.p. function itself. We will denote the one in ([T4| by p". It corresponds to what we 
will refer to as "the augmented system" , whose units are of two kinds: the original 
system's components forming a set [n] and the set of binding factors denoted by 
[rn]. Now equation (fT4|) becomes 

r5=p'^(T°,...,T°,Qi,...,Q™). 

Under a natural assumption that intrinsic components' lifetimes are indepen- 
dent of the bounding factors, the augmented system falls into the class of Bayes- 
dependent systems where the bounding variables play the role of dependence- 
inducing factors. To apply the relevant formulas of Sections 2 and 3 to this case 
the following notation will be needed. 

• Joint reliability function of bounding variables: 

Rb{t) = Pr(Qj > Vj e H); 

• Component subsets of the augmented system: ^ © B, A C [n], _B C [m]; 

• v"". the "u" -function corresponding to p"^ . 

Now formulas (HI]) and yield that 



(15) Rs{t) = I] ™«»(^©5)-R6(e^*) H-^^W' 

A<Z[n]B<Z[m] iGA 

poo 

(16) MTTFs - Y 51 ™"''(^®^) / Me°B)Y[Mt)dt, 

AC[n] BC[m] "^^ i£A 

A practically important special case arises under additional assumptions that 
the original components' intrinsic lifetimes have exponential distributions and that 
bounding factors are independent of each other. Then 

Rb{e°/) = n ^^■(^) 

jeB 

and formulas ([T5|) and (fTB]) yield the following: 

= rn„{A®B)e-^^'l[R,{t), 

ACln] BC[m] jeB 

/"OO 

MTTFs = 51 "'^(^®-^) / e-^^' Y[Rj{t)dt, 

AC[n] BC[m] ° j£B 
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where, as before, Xa = XlieA '^i- 
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